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1 Introduction 

The Bose-Einstein correlations in multiple particle production processes were 
discovered almost fifty years ago [87 , [88J. Since then many old results got 
forgotten. For instance now most people believe that what is being observed 
in multiple particle production processes is the Hanbury Brown and Twiss 
effect, in spite of the very clear proof to the contrary given by Kopylov 
and Podgoretsky |114j (cf. |119j for a more recent brief discussion of this 
point). Concepts are renamed. E.g. what is now called Bose-Einstein or 
HBT correlations used to be the GGLP effect. Results are rediscovered and 
priorities are ascribed at random. It seems, therefore, useful to summarize 
and review critically the old results. We chose the first thirty years, which 
corresponds to little more than two hundred papers. Even so, this is a vast 
subject and I will be grateful for comments and remarks, which could help to 
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improve the text. After 1989 the explosion in the number of publications on 
the subject has been such that now the number of papers can be estimated 
to be over one thousand. Around 1990 the first big review papers devoted to 
Bose-Einstein correlations appeared |134j . |123j . |4*£] . For obvious reasons, 
however, they concentrated on the problems which the authors considered 
topical at the time of writing and did not have the hindsight we have now. 

This review is chronological, but not historical. A historical review would 
have to include all the wrong and/or sterile ideas which influenced the re- 
search workers in the field. No effort in this direction has been made. 

When reviewing the work of many people, done over a long period of time, 
the choice of notation becomes a problem. We have chosen a simple, uniform 
notation without trying to keep track of the notations of the various authors. 
Our notation will be introduced when necessary, but here we summarize the 
main points. The four-momentum of a particle will be denoted p or p' and 
its positions in space-time x or x'. If the particle is on its mass shell, the 
time component of p is E p = a/p 2 + m 2 . For a pair of particles we introduce 



K = ^(pi+p 2 ), K' = -{p[ + p' 2 ), q = pi-p 2 , q' = p[-p' 2 , U: 



Q 
X 



-Q 2 



(2) 



(xi + x 2 ), X' = -(x[ + x' 2 ), Y = xi - x 2 , Y' = x[ - x' 2 ,(3) 



with the obvious notation wQ 2 = Q. For single particles 



K i = g(Pi +Pl)> Q f i=Pi-Pi> x i = Y i = x i- X 'v ( 4 ) 

The subscript is changed from 1 to 2, if the position in space-time and the 
four-momentum of the particle are x 2 ,P2- Sometimes, when only one particle 
is considered and there is no risk of confusion with the previous formulae, 
the subscripts will be dropped. 

Let us note for further reference two useful identities 

PiXi+p 2 x 2 = 2KX-^qY, (5) 
P1X1 - p 2 x 2 = KY + qX. (6) 

and their analogues for the primed and subscripted parameters as well as for 
the space vectors. 
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We will use pseudo density matrices p(p 1; . . . , p n ; p[, . . . , p' n ). Their diag- 
onal elements yield the momentum distributions, but the relations are model 
dependent and in general more complicated than for the standard density 
matrices. The pseudo density matrices are hardly ever normalized to unity. 
Moreover, in the GGLP model discussed in the next chapter their diagonal 
elements have to be integrated over phase space, with energy momentum con- 
servation imposed, in order to give the momentum distributions (see text). 

The radius R of the interaction region means different things in different 
papers. For a Gaussian distribution of sources 

\J (r 2 ) = ^3/2aR, but both a = 1 and a = 1/2 are being used. For sources 
distributed on the surface of a sphere (r 2 ) = R, while if the sources are 
distributed uniformly over this sphere y/ (r 2 ) = a/3/5-R. One should also 
keep in mind that the average of the square of the distance between two 
points is 2(r 2 ), while e.g. a/ (x 2 ) = y/l/%\J (r 2 ). These ambiguities have 
been discussed and illustrated by many examples from the literature in ref. 

In the next section we review the two papers which created the field. 
Then we present the subsequent results divided in ten years periods. Finally 
in the last section we summarize the highlights of this development. 

2 Beginnings 

2.1 Discovery 

Bose-Einstein correlations in multiple particle production were discovered ac- 
cidentally jHH] as a byproduct of an unsuccessful attempt to find the p meson, 
which had just been predicted by Frazer and Fulco [82J. The group of the 
Goldhabers at Berkeley was studying annihilations of antiprotons of momen- 
tum 1.05 GeV/c in a propane bubble chamber. Since the predicted p meson 
had isospin one, they hoped to find a bump in the invariant mass distribu- 
tion of the unlike-sign pairs (tt + 7t~), and no bump in the mass distributions 
of the like-sign pairs (n ir ). They analyzed the 2ir + 2ir~ and 37r + 37r~ semi- 
inclusive channels 1 [B2J , but instead of finding the expected bumps discovered 
that the distributions of the angle between pion momenta, as measured in 
the pp cms frame, strongly depended on the pion charges. For like-sign pion 

1 i.e. channels with an unspecified number of tt°-s 
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pairs the distributions in cos 6^ were almost flat, while for unlike-sign pairs 
large opening angles were more probable as expected from momentum con- 
servation,. In fact, the distributions for the unlike-sign pions were even more 
peaked than required by phase space, so that the combined distributions of 
cos8 n7T for all the pion pairs at given multiplicity agreed within errors with 
calculations from the statistical model with Lorentz invariant phase space. 
In order to characterize the distributions of the opening angle, the group 
introduced for each distribution a parameter which was being very popular 
for many years, 



where n > (n<) is the number of pion pairs with the c.m.s. opening angle 
greater (smaller) than 90°. Thus, the experimental result was that r ) l%ke < 
^unhke ^th the value of 7 calculated from Lorentz invariant phase space 
falling in between. 



It took over a month and the help of an eminent theorist to understand that 
what was being observed were correlations due to Bose-Einstein statistics 
jSHJ- The resulting paper written by the Goldhabers, Lee and Pais [BB] (fur- 
ther quoted GGLP) became a standard reference. According to the SPIRES 
data base, which omits citations previous to year 1973, this paper has been 
quoted over 300 times. For many years the effect was being referred to as the 
GGLP effect. The realization that an analogous effect had been discovered 
earlier in astronomy by Hanbury Brown and Twiss jUH] and applied to the 
determination of the angular radii of stars came only much later [94], |%9*] . 

The starting point for GGLP was the statistical model with Lorentz in- 
variant phase space |145j . According to this model the distribution of the 
particle momenta for a final state produced in a pp annihilation, containing 
n pions and nothing else, is given by the formula 



Here a is a constant, P is the initial total momentum four- vector and the 
momentum four- vectors of the final particles are (E Pi , pj, i — 1, . . . , n. The 
factor C contains a normalization constant and an isospin factor dependent 
on the numbers of positive, neutral and negative pions (n+, no, nJ). It will 
not be discussed here any further. Today it is well known that prescription 



n < 




2.2 The GGLP paper 




(9) 
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Q is very unrealistic, but at the time of GGLP, when beam energies were 
low and for comparison with experiment the formula was integrated over all 
the degrees of freedom except at most one, agreement with data was usually 
reasonably good. 

The point made by GGLP was that formula (jHJ) contains no correlations 
between the electric charges of the pions and their momenta. Consequently, 
the distribution of opening angles for pion pairs, as predicted by the model, 
does not depend on the charge of the pair, which is in violent contradiction 
with the data |HZ|. They proposed to replace the factor l/n\ by a product of 
three factors, one for pions of each charge. In view of future generalizations 
we will present their argument using operators and matrices instead of the 
wave functions they used. By analogy with the formulae for the density 
matrices in the momentum representation, we are looking for a replacement 

1/n! ->IJp fc (p 1 ,...,p nh ;p 1 ,...,p n J, (10) 

k 

where the product is over the pion charges k = +, 0, — and the tildes are 
there to remind that pk do not have all the properties of textbook density 
matrices. We will call the matrices p pseudo density matrices. This includes 
a variety of matrices similar to, but not identical with, the density matrices. 
Let us define a single particle density operator 



p = J (fx|x)p(x)(x| (11) 
with the usual normalization condition 

d 3 xp(x) = 1. (12) 

Function p(x) must be real, because operator p is hermitian. This single 
particle state is an incoherent superposition of pure states, each correspond- 
ing to a given production point x. GGLP considered two weight functions 
jo(x). Following the statistical model they put p(x) =const within a spherical 
volume of radius Rs and zero outside. They also introduced a "Gaussian- 
shaped volume" meaning that p(x) ~ exp[— x 2 /2A]. Comparing the matrix 
elements p(Pi) P2> Pi> P2) calculated within their model (see further) for the 
two choices of p(x) they found that the results were similar, with the closest 
agreement for 

#5 = 2.15^. (13) 
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Therefore, they chose the Gaussian as it leads to simpler calculations. The 
density matrix in momentum representation, which corresponds to ([lip, is 

P(P; PO = (PlPlpO ~ J d 3 xp(x)e-^ = (e^) . (14) 

The corresponding momentum distribution can be obtained by putting q = 
and is constant, as was to be expected from the uncertainty principle. Let 
us note that p(p; p') is a single particle density matrix, but not necessarily 
the one describing the system being considered. 

Let us define an unsymmetrized n-particle density matrix as a product 
of the single particle density matrices ()14j) 

n 

P U {PI, • • • , Pn> P'l> • • • > Pn) = II p ( p i> ( 15 ) 

i=l 

Taking the diagonal elements and imposing energy momentum conservation 
one could "derive" in this way the statistical model ®. GGLP, however, 
proposed to symmetrize for each charge of pions. For k = +, 0, — their 
assumption, written in a more general form |l()5j , is 

Pfc(Pl) • • • ! Pn fc > Pi; • • • 5 Pn fc ) = ~~~| Pk (Ppi) • • • > PPn fc ) PQ1) • • • j PQn fc )> 

fc ' P,Q 

(16) 

where the summation extends over the n^! permutations P of the momenta 
Pj and over the n^! permutations Q of the momenta p^. This formula can be 
rewritten as 

Pfc(Pi; ■ ■ ■ j Pn fc > Pi) • • • ) Prife) = / ] Pfc (Pi) ■ ■ ■ ) Prife) PQ1) ' ' ' ) PQn fc 

Q 

GGLP considered no more than two pions of the same sign, because they 
assumed in the calculations that every event with more than one ir° contains 
exactly two 7r°-s. Moreover, they needed only the distributions of momenta 
i.e. the diagonal elements of the pseudo density matrices. Thus, their recipe 
was as follows. For nf. — 0, 1 there is no dependence of the factor on the 
momenta of the pions. For = 2 

Pfc(Pi,p 2 ;Pi,p 2 ) = p(pi;pi)p(p 2 ;p 2 ) + Ip(pi;p 2 )| 2 , ( 18 ) 
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where the hermiticity of the single particle pseudo density matrix has been 
used. Substituting expression (fHj) one finds: 



(19) 



Pfc(Pi,p 2 ;Pi>P2) = 1 

which for a Gaussian p(x) gives 

P fe (Pi,P 2 ;Pi,P 2 ) = l + e- A(Pl - p2)2 . (20) 

As seen from the derivation, the momenta in the exponent are three- moment a. 
GGLP, however, replaced them by four-momenta, in order to simplify fur- 
ther the calculations. This makes the formula Lorentz invariant and does 
not change much their numerical results. Incidentally, it improved agree- 
ment with experiment. Let us follow GGLP and introduce their "correlation 
function" 

tfj( Pl ,p 2 ) = l + e- XQ \ (21) 

where Q 2 = — (pi — £> 2 ) 2 - This formula has later become very popular. Let 
us stress, however, that in the GGLP model it was a factor in the integrand 
of the phase space integral, while later it was being used without further 
integrations. For the channels 27T+27T", 27r + 27r~7r° and 27r + 27r"27r° the GGLP 
proposal amounts to replacing in formula (jH} the factor 1/n! by 

^(Pl ) P 2 )V'(P3,P4)V'(P5 ) P6), (22) 

where pi and p 2 are the four-momenta of the positive pions, p% and p^ are 
the four-momenta of the negative pions and 4>{P5iP§) = 1 unless uq = 2. For 
n = 2: ^{pbiPo) = ip(p5,Pe,), where ps and p§ are the four-momenta of the 
two neutral pions. The results, after averaging with suitable weights over 
n , were found to be in good qualitative agreement with the corresponding 
experimental results from |87j . Also the radius of the interaction region Rs 
came out between 0.7 fm and 1 fm which is the expected order of magnitude. 
Quantitatively, however, ^ unhke was predicted too small and r y hke too large 
whatever the choice of Rs- Thus the model underestimated the difference 
between the like- and unlike-sign pion pairs. 

Besides their calculations GGLP made some important general remarks. 
They realized that "an adequate model should at the same time give a rea- 
sonable account of all combined aspects of the annihilation process" and that, 
therefore, one should look for other evidence for or against the model. They 
predicted that for any given exclusive channel the effects of Bose-Einstein 
statistics should decrease with increasing energy. They worried about the 
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multiplicity distributions, because the volume |vri?|, which they found, was 
significantly smaller than the volume which in the statistical model governs 
the multiplicity distribution. They stressed that the study of Bose-Einstein 
correlations may give valuable information about the reaction mechanism, 
but that "results of this study should not be construed to imply that dynam- 
ical effects (such as, for example, n-n interactions) are definitely negligible" . 
In fact they suspected that the "it isobars" , or in today's language resonances, 
may be important. 

3 First decade (1961 — 1970) 

In the sixties the study of Bose-Einstein correlations was gradually gaining 
popularity. The conceptual framework was still that of the GGLP paper, 
though an isolated attempt to study three- and four-body correlations [ZOj 
should be noted. Usually, people plotted the distributions of cos9 n7T , calcu- 
lated the parameters 7, or some other parameters which could be expressed 
in terms of the parameters 7, and supported, or very rarely criticized, the 
implications of the GGLP paper. The effect found in jHZj was confirmed for 
other pp interactions. A compilation from 1967 jOZ lists 7", 7' pairs for 25 
different (exclusive) channels and/or energies. The effect was found in ir + p 
7i~p |lU3j . pp jnni and K + p jSH] interactions. The effective attraction in 
momentum space was found for 7r°7r° pairs [TS], but not for K + tt ± pairs [6*H] . 
The results for the ir°-s and the K + tt + pairs eliminated the possibility |140| 
that the effect for 7r ± 7r ± pairs was due to the fact that these pairs are "exotic" 
i.e. unable to form resonances, and not to the fact that they consist of iden- 
tical bosons. A related possibility was that the effect is due to resonances. 
Here the results were not so clear (cf. eg. the review [150J and references 
given there), however, cutting out the resonances did not change much the 
measured 7-s [HE]- Moreover, a particularly strong effect was observed for 
7i + p interactions at 8 GeV/c with nine charged pions and a proton or ten 
charged pions and a neutron in the final state, where the energy per particle 
is too low for a significant production of resonances • Thus Bose-Einstein 
correlations remained as the only plausible explanation. 

Some qualitative implications of the GGLP model got confirmed. A com- 
pilation of the data for annihilations showed [67J that with increasing energy 
the effect, as measured by the parameter 

° = ( r + i)(y + i)' (23) 
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decreases rapidly. A rough fit was [HZ] 



C 



73.9 




(24) 



<- 64 (n + 1) 



1.54 



where n c and hq denote the numbers of charged and neutral pions respec- 
tively. The centre of mass energy ^/s is here expressed in GeV. The un- 
certainty of the normalizing factor 73.9 was about fifty percent, while the 
uncertainties of the other three parameters were at the ten per cent level. 
For events with many particles in the final state it was possible to show that 
the effect occurs in every event and not only in some special ones The 
effect is enhanced when one chooses pairs of pions with similar values of 
momenta |pj i.e. with similar energies ITfj . 

Quantitatively, however, the model did not work. The conclusion in joT] 
was that "The model of Goldhaber et al. is in violent disagreement with the 
data". The main argument was that at low energy per particle the values 
predicted for C are much too small. Moreover, assuming a constant inter- 
action radius, R — 1.3 fm, the authors found that for the exclusive process 
pp — > 3iy + 3iy~ the model grossly underestimates the energy dependence of 
parameter C. Incidentally, [HZ] seems to have been the first calculation with 
three identical particles in the final state. The statistical model contains no 
peripherality and consequently with increasing energy per particle its pre- 
dictions for the single particle distributions become completely wrong. The 
GGLP model did nothing to correct that, but the point was that also its 
predictions for the correlations were wrong. One tried to improve the start- 
ing point by replacing the statistical model by the uncorrelated jet model, 
which includes peripherality, or by the CLA model, which included moreover 
some multiperipheral correlations, but this did not help [68J. The situation 
was summarized at the 1-st multiparticle conference in Paris |15Uj . To put it 
short: the GGLP model was good enough to convince people that the effect 
is due to Bose-Einstein statistics, but not good enough to make quantitative 
predictions. 

4 Progress in the seventies 

4.1 The approaches of Kopylov and Podgoretskii, and 
of Shuryak 

In 1971 an important series of papers by Kopylov, Podgoretskii and collab- 
orators (further quoted KP) begun to appear jM]. The approach of KP 
differed significantly from that of GGLP. There was a major difference in 
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motivation. The purpose of GGLP was to explain the charge dependence 
of the distributions of the opening angles for pion pairs. They noticed that 
the radius of the interaction region, obtained as a byproduct, was of the 
expected order of magnitude, though probably too small, but this remark 
was not followed. According to a detailed review [48J only two experimental 
papers previous to 1976 gave values for the radius of the interaction region 
(jHZj and jZHj). Looking up these papers, one sees that in both cases the radii 
were given only to convince the reader that the authors worked hard before 
recognizing that they are unable to fit the data. KP were familiar with the 
work of Hanbury Brown and Twiss [OH] and their purpose was to use quan- 
tum interference, for identical fermions as well as for identical bosons, to 
get information about the lifetimes and extensions of various sources. They 
considered a variety of sources |112j : resonances of particles and of nuclei, 
highly excited nuclei where both lifetimes and shapes were of interest and 
interaction regions in multiple particle production processes. KP did not 
compare their results with experimental data. 

The GGLP approach did not include time. It corresponds to a picture 
where all the pions are emitted instantaneously and simultaneously. Then, 
when it exactly happened is irrelevant. The sources of KP were time de- 
pendent and extended in time. They considered it to be the main difference 
between the two approaches [03] an d claimed to have "a more correct theoret- 
ical approach" |112j . Another difference was technical, but of great practical 
importance. GGLP had reasonably simple formulae for the integrands of the 
phase space integrals. The integration over phase space, however, remained 
to be done. In practice this limited their calculations of the physical distribu- 
tions to exclusive channels with no more than six particles. KP assumed that 
the constraints of energy and momentum conservation have little effect on 
the distributions of the momentum difference (q , q) = (E P1 — E p2 , p 1 — p 2 )in 
the region of small |q|. Thus, their results could be directly compared with 
the measured cross-sections. This eliminated the necessity of evaluating the 
phase space integrals and made the predictions applicable also to inclusive 
processes. Obviously, it was a reasonable assumption only for sufficiently 
high multiplicities for exclusive processes, or for sufficiently high energies for 
inclusive processes. The KP results for multiple particle production were 
sketched in |lllj and then described in more detail in |112j . In ref. |113j 
the theory was reformulated using a more powerful formalism borrowed from 
optics and some generalizations were given. A short summary, including im- 
portant new ideas, was given in |lU9j . Reference |114j explained the relation 
between the work of GGLP and that of Hanbury Brown and Twiss. 

Both GGLP and KP considered an incoherent superposition of pure pion 
states produced by sources excited at some common time, say t — 0, at vari- 
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ous points of space x and emitting one pion per source. The GGLP sources, 
however, emit the pions instantly and without changing their locations. The 
KP sources are much more general. The amplitude of a pion emitted by a 
source was assumed to be the solution of the Klein-Gordon equation 

d 2 \ 

V 2 - — - m 2 ) A(x) = -4ttJ(x), (25) 



dt 2 

where m is the pion mass, x is the four-vector x, t and J(x) is the (classical) 
source. Further there will be many sources contributing incoherently, but let 
us start with just one. The equation simplifies when the Fourier transforms 
in time: 

/OO /' X, 

dtA(x)e i£t ; J(x;e)= / dtJ(x)e iet ; (26) 
-OO J — OO 

are introduced. Note that, according to the general rules of quantum me- 
chanics, e can be interpreted as the energy of the pion. The transforms 
satisfy the equation 

(V 2 -p 2 )A(x;£) = -4vrJ(x;£), (27) 

where p 2 = e 2 — m 2 and, therefore, it is the square of the momentum of the 
pion. The solution of this equation is well known and described in detail in 
|113j . but we will need here only a special case. Let us assume that all the 
source functions J(x) are negligible outside a small region of space V and 
that the detector is at a large distance ro from the centre of region V. Large 
means here much larger than the diameter of region V. Then, in the vicinity 
of the detector, the (approximate) solution is 

A(x;e) =A(p)^—^- A(p) = [ d 4 xJ(x)e ipx - (28) 
r D J 

where p is the pion energy-momentum four- vector (p, E p ) with the momen- 
tum p oriented from the centre of region V to the detector and |p| = a/p 2 . 
Note that all the dependence on the source is contained in the factor A{p). 
The other factors are irrelevant from the point of view of symmetrization and 
will be further omitted. Thus in general one source, say the one labelled i, 
contributes to the single particle pseudo density matrix a term Ai(p)A*(p'). 
In general this product depends on some parameters and has to be averaged 
over them. The result should be summed over all the sources. Note that the 
contributions from all the space-time points where J(x) is significantly dif- 
ferent from zero are summed coherently (adding amplitudes) in the formula 
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above, while the averaging and the summation over sources correspond to 
incoherent (adding probabilities) superpositions. Finally 

p(p,p') = J2(Mp)aKp'))- (29) 

i 

Formulae of this kind, which are more general than the KP model, were 
introduced by Shuryak |141[ 1142] . Shuryak's results can be derived starting 
from the density operator 

p{t)= f dto f <e-^- t °)|J(to))(J(t / )|e + ^- t o). (30) 



Actually it is a pseudo density operator, because its trace does not have 
to be equal one. Further we will take t — > oo. When J(x) = (x|J(i)) is a 
classical current this leads to the KP theory. J(x), however, can be any single 
particle production amplitude (cf. e.g. |139j ). This formula, as it stands, is 
for one coherent source. In order to make a realistic model it is necessary 
to sum over the sources and, sometimes, to average over some parameters 
of the sources. Including that, introducing, after the first and before the 
last factor of the integrand, unit operators built from the eigenstates of the 
position operator and converting to a pseudo density matrix element in the 
momentum representation we find 

(31) 

Evaluating the matrix elements of the time-evolution operators and omitting 
the uninteresting factor exp[it(E p > — E p )] we get 

m-J*I«($ww))<~<>, m 

or equivalently in other variables 



(33) 

This formula is the starting point for much modern work. 
KP proposed 

J[x) = r- 1/2 5(x - x - vt)e- tEt -^6(t). (34) 
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This corresponds to a classical point source created at time t — at point 
x = x and moving from that point on with constant velocity v. The source 
has energy E and decays emitting a pion according to an exponential decay 
law with average life time r. Thus the parameters xo, v, E, r depend on the 
source label i. The integrations necessary to get A(p) are in this case trivial 
and one finds using (j2H|) 



-iq-xo 



PM) - E ^ (i + lT{E _ e + p . v)) (i _ lT{E _ £ , + p > . v)) ) • ( 35 ) 

Here KP assumed that the averaging should be done over the source energy 
E and that this can be reduced to an integration over E from minus to 
plus infinity. The integral is easily done using the method of residues. The 
sources are labelled by the points x of their creation. Thus, the summation 
over % reduces to an integration over space with some weight p(x ). Finally, 
omitting an irrelevant factor 27r, 

/d 3 x p(x )e-^ x » 
1 - tr(q - q ■ v) 

Usually one puts v = [109] so that 



|P(P'P)I = 1 , _2^2 • ( 37 ) 



1 + r 2 g 2 

The symmetrization is done as in GGLP. Thus for the two-body correlations 
formula (jl8j) remains valid. Since p(p, p) = 1, the first (non interference) 
term in (j!8|) remains equal one. The second (interference) term is modified 
by the factor (1 + r 2 ^^)" 1 . From the experimental point of view this is a 
beautiful formula. From the dependence of the correlation function on the 
vector q at fixed q one can determine the space distribution of the sources 
just like in the GGLP model. Studying the dependence on q at fixed q, 
however, one can additionally measure the life time of the source. 

An obvious objection to formula (}36|) is that a density matrix depend- 
ing on the difference of the momenta only is grossly unrealistic. GGLP 
did not have this difficulty, because they kept only the states allowed by 
energy-momentum conservation and the necessary projection introduces a 
dependence on pi + p 2 . KP [113] proposed the following way out. Suppose 
that the incoherent sources are not point-like in space, but smeared. Thus 

d 4 x s e- tH ^-^\^ s }p(x s )^ s \e +lHo ^\ (38) 
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Here the four- vector x s labels the sources. It gives the space and time position 
of the source at the moment when the source got created, or of any other 
point in space-time which defines unambiguously the source. Function p{x s ) 
gives the distribution in space-time of the incoherent sources, thus, for any 
function f(x s ) the average over the sources is 

(f(x 9 )) = J d 4 x s p(x s )f(x s ). (39) 

Hq is the free particle Hamiltonian. The exponentials give the time depen- 
dence of the state vectors \ip s ) and their dependence on the time component 
of x s . The state vectors \ip s ) are equivalent in the sense that 

<x|^)=#K-x a ). (40) 

Thus all the incoherent sources are related by rigid shifts in space-time. It 
is convenient to use a picture where the density matrix does not depend on 
the time t. This can be achieved either by using the time dependent momen- 
tum eigenstates e~ lHot \p), or equivalently by using the interaction picture, 
where each state vector gets an additional factor e +lHot , which cancels the 
t-dependent factors both in the density operator and in the time dependent 
momentum eigenstates. Let us define the function 

A(p) = J rf 3 xe- ipx ^(x). (41) 

Then, taking the matrix element of the density operator and dropping some 
constant factors, 

p(p 1 ,p 2 ) = A( Pl K(p 2 )<e^>. (42) 

This formula has several interesting features. The diagonal elements, which 
yield the single particle momentum distribution are 

p(p,p) = |A(p)| 2 . (43) 

Thus, any single particle distribution can be reproduced perfectly. Function 
A(p) is related to the properties of a single source defined by the state vector 
but tells us nothing about the distribution of sources in space-time. On 
the other hand, for the correlation functions defined in the next section the 
factors A(p) cancel and the formulae are as if the sources were point-like. 
Therefore, the correlation function depends only on q and not on the vectors 
Pi,P2 separately. This is similar to the situation in the GGLP model. Note, 
however, the change in the interpretation of the dependence of function p(x s ) 
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on the space vector x s . Now it gives the distribution of the positions of the 
sources. Since the sources are smeared in space, the actual region, where the 
pions are produced is larger than would follow, if the function p(x s ) were 
interpreted as the GGLP function p(x) was. What is more, the function of 
four variables p(x s ) cannot be determined unambiguously from measurements 
of the three-momenta. Thus unavoidable model dependence comes in. 

Another less important, but much publicized, difference between KP and 
GGLP was in the assumption about the space distribution of sources. KP 
pointed out that the interaction region is likely to be opaque and that, con- 
sequently, only the sources on its surface can contribute. They assumed that 
the sources are created on the surface of a sphere of radius R, but that their 
radiation satisfies Lambert's law taken to mean that this surface can be re- 
placed by a circular disc of radius R perpendicular to the direction from the 
interaction region to the detector 2 . Thus for the integral in they got 



Here is the length of the component of q in the plane of the disc and 
Ji is a Bessel function. The integral is most simply done by expanding 
the exponent in the integrand and integrating term by term. Cocconi [HS] 
proposed to replace this expression by a Gaussian, according to 



Actually, for not too large values of the argument, the right hand side is a 
very good approximation to the more complicated expression on the left-hand 
side. This short reference does not give justice to the importance of reference 
[63J. Cocconi derived independently many of the KP results. He proposed 
alternative interpretations. E.g. in his picture cr was the thickness of the 
"photosphere" i.e. of the spherical crust from which the pions were emitted. 
Last not least, he contributed much in discussions with experimentalists, so 
that ref. jHS] has got more citations than any of the KP papers and even 
sometimes the whole approach was referred to as the Kopylov, Podgoretskii 
Cocconi model. 

2 This is a hand-waving argument. An exact evaluation of the corresponding integral 
can be found in ref. |115j . It yields the result 144|l in the limit ct ^> R. This condition 
was clearly stated in refs |112| . Unfortunately, it is hardly ever realized in practice. 

Another case when the formula is valid is for q\\ = 0, where q\\ is the component of vector 
q along the direction of the total momentum of the pair. 




(44) 




(45) 
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For the distribution in the difference in momenta KP got |lU9j 



(Pqxdqo qrR 
Kopylov |l()9j rewrote this in the form 



d 3 a \2Ji(q T R) 
; — ~ 1 + — 




(46) 



d 3 a 



d 3 a 



) 



I(q T R) 
1 + g V 



) 



(47) 



d 2 qxdqo 



d 2 qxdqo 



1 + 



off 



Note that this formula would not follow, if the second factor on the right-hand 
side depended on p 1 + p 2 , because the two-particle pseudo density matrix 
yields d 6 a/dp 6 and after integration over P! + p 2 factorization would not hold 
any more. Function I(qj-R) can be obtained by comparison with the previous 
formula, or from some other model. The important distinction is between the 
cross-section on the left-hand side, which is the physical one containing all 
the interference effects, and that on the right-hand-side, which is calculated 
with the interference effects turned off. In the GGLP and KP models this 
would be up to a constant factor p(p 1 , p 1 )p(p 2 , p 2 ), but the formula is more 
general than that. Then the "off' cross-section, which was soon renamed 
by experimentalists "background", should be taken from experiment. The 
question is how? GGLP used for comparison with like-sign pion pairs the 
unlike-sign pion pairs. These have certainly no Bose-Einstein correlations, 
but they have other unwanted correlations due e.g. to the formation of 
resonances in the tt + it~ system. Kopylov |109| proposed to take like sign 
pairs, but with each pion taken from a different event 3 . This background 
got later the name "mixed" and became quite popular, though it eliminates 
too many correlations e.g. the correlations due to energy and momentum 
conservation or to the production of jets. 

The relation between the work of GGLP and that of Hanbury Brown and 
Twiss was clarified in ref. |114j . The authors developed a more general model. 
It introduces two scales of time, or energy: the life time of the emitting 
source, or equivalently the width of the energy distribution of the produced 
particles, r ~ ^ and the resolution of the detector in time, or equivalently in 
energy, At ~ -. In astronomy the life time of the excited atomic state, which 
emits the radiation, is much shorter than the time resolution of the detector. 
Therefore, one can neglect r compared to At. In particle physics, on the 
other hand, the width in energy of the source is much larger than the energy 
resolution of the detector. Therefore, the opposite limit is justified: one 
can neglect 7 as compared to T. These two limits correspond respectively 
to the models of Hanbury Brown and Twiss and to that of GGLP. They 

3 According to ref. |11U| this was an idea of Podgoretskii. 
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are completely different limits. E.g. the HBT results can be derived from 
classical physics without ever mentioning quantum physics. Nevertheless, 
the acronym GGLP got in the seventies gradually replaced by the acronym 
HBT. The analysis from |114j implies that the motivation must have been 
sociological rather than physical. In this review we use the neutral BEC, 
which is also quite popular. 



4.2 Correlation functions 

Since the study of various correlations was a popular subject in the seventies, 
it soon became common to express formulae like (|47J) in terms of two-particle 
correlation functions, cf. e.g. |H], defined by 



Here 

1 d 3 a , , 1 d 6 a 

= = ;*55r (49) 

are respectively the single particle and the two-particle distributions of the 
particles studied, a is the integrated single particle cross-section, and d 3 p 
means either dp x dp y dp z or dp x dp y dp z / E p . Af is a normalization constant. 
Theorists usually either put it equal one, or interpret it as the ratio of the 
normalization of the denominator to the normalization of the numerator. In 
the latter case, for inclusive processes they put Af = j^zfjj- The resulting 



confusion has been reviewed in ref. [128J. Experimentalists more often use 
this factor to make the correlation function tend to one with increasing Q 2 . 

The correlation function is a nice, well-defined object, but for extracting 
information about the production region it is often "extremely impractical" 
as put later by experimentalists [148] . Consider for example a process, where 
the final state consists of two narrow jets back to back. Let the orienta- 
tion of the axis of this pencil-like structure have an isotropic distribution. 
Then function N(p 1 ,p 2 ) nas strong maxima for 8 n7T = and 6 n7T = tt, while 
function N(p) is spherically symmetric. The correlation function also has 
a strong maximum at 6 nn = (and another one at 6 n7T = tt), but this tells 
us little about Bose-Einstein correlations. On the other hand, if for each 
event the coordinate system is rotated so that the axis of the jets is along 
the z axis, also function N(p) acquires maxima at 9 n7T = and 9 n7T = tt. 
The maxima in the denominator of (}4*8j) largely cancel the dynamical part of 
the maxima in the numerator and Bose-Einstein correlations can be studied, 
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though usually some additional selections to improve the analysis are intro- 
duced ( c.f. e.g. Other examples, where the correlation function (pIHj) is 
in practice useless for gaining information about the interaction region, have 
been described in ref. |95j . Formally, the way out is to define the correlation 
function Co(j»i,]92) for the background, i.e. for the fictitious case considered 
by KP: the same process, but with the Bose-Einstein correlations switched 
off. Then the reduced correlation function 

C(Pi, P 2 )/Co(Pi, p 2 ) = 1 + R(Pi, P 2 )- (50) 

is a good replacement for the correlation function ()48j) . In the following we 
will mostly discuss the function i?(p 1 ,p 2 ) defined by this formula. It corre- 
sponds to the interference term in formula ()47|) and should tend to zero, when 
Q 2 becomes so large that no more Bose-Einstein correlations are expected. 
When the background two-particle distribution is a product of the standard 
inclusive single particle distributions, Co(p 1 ,p 2 ) = 1 and 1 + R(p 1 ,p 2 ) re- 
duces to the correlation function C(p 1 ,p 2 ). This is the case for the mixed 
background, while usually it is not the case for the background from unlike 
sign pion pairs. It is also not the case when the inclusive single particle dis- 
tribution is modified by some selection or transformation like in the two-jet 
model mentioned above. Since it is usual to call 1 + R(p l ,p 2 ) a correlation 
function, one could argue as follows. The uncorrected experimental corre- 
lation function C uncorr (p 1 ,p 2 ) should be corrected for the unwanted corre- 
lations by dividing it by Co(p 1; p 2 ). Then 1 + _R(p 1; p 2 ) is this corrected 
correlation function. In the seventies theorists usually put Co = 1 and ex- 
perimentalists chose Co constant and such as to make R{pi, p 2 ) tend to zero 
with increasing Q 2 . Later the sophistication was greatly increased, but the 
best way of choosing Co(p 1; p 2 ) is even today controversial. 

4.3 Further experimental results 

On the experimental side some more qualitative implications of the GGLP 
picture got confirmed. The effect was found in nucleus-nucleus collisions jHH| 
and for K®K® pairs p4"] . It was shown [77] that within the experimental 
errors the GGLP effect could be interpreted as a reflection of the difference 
in the distributions of the masses of pion pairs for pairs of like-sign and 
unlike-sign pions. Actually GGLP proposed that the weight factor, which 
equals one for unlike-sign pairs, should depend on Q 2 = —{p\ — P2) 2 for like- 
sign pairs. Since, however, Q 2 = m 2 n — 4:fn 2 , the two statements are equiv- 
alent. Three- and four-body correlations were studied jUJ, |118j . They 
were found strong, but at that time could be fully explained as reflections of 
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the GGLP two-body correlations. Gradually the groups begun to determine 
the radii and life times of the interaction regions as suggested by KP. Re- 
sults were obtained for K + p scattering [HI], [021, f° r n P scattering [22], 
jHHj , [ZD] , for pp [IB] , [H] , [M] , |122j and pn [HOI scattering , for pp scatter- 
ing jH], [Hj, [ZH|, for p-nucleus (Be,Ti, W) [H| and tt~ - C [23! scattering 
as well as for nucleus-nucleus scattering jHH] {Ar + Bal 2 and Ar + P6 3 4 ), 
(C + Ta). The background distribution was usually identified with the 
distribution for unlike sign pairs. Only for scattering on nuclei the mixed 
background proposed by Kopylov was found more convenient, because using 
it one could use the negatively charged particles only and thus, avoid the 
contamination by protons. Sometimes results for two different backgrounds 
were compared in order to show that the uncertainty about the background 
distribution does not affect too much the results. The normalization of the 
background distribution with respect to the distribution of like-sign pairs was 
chosen so that the numbers of pairs in both samples coincided, or so that 
the distributions became identical for large values of Q 2 . Note that these 
two normalizations are not equivalent. There were various other sources of 
uncertainty. Statistics was poor and the results depended on the binning 
chosen. The fits, as judged by the \ 2 test, were often poor. The maps of 
X 2 (-R, r) exhibited large valleys with secondary minima. Formula (|37J) is not 
covariant. Some groups used it in the overall cms, others preferred the cms 
of the charged pions. On that were superposed the ordinary experimental 
uncertainties e.g. particle misidentification and limited momentum resolu- 
tion. The theory was not very robust either. E.g. when the authors of ref 
|122j studying pp annihilations got i? ~ cr « 3fm and, not very surprisingly, 
considered this result too large, they replaced the assumption that v = 
by the assumption that there are two sources moving with velocities v and 
— v respectively, added some phenomenological assumptions about v and got 
from the same data cr « 0.7fm and R ~ 0.7 — 0.8fm. Another degree of free- 
dom, which was not exploited by experimentalists however, was to assume 
that there is a scatter of the time when the sources got excited. This gives 
another factor dependent on q |112j . 

In spite of these uncertainties the results were not unreasonable. Results 
from K®K® pairs [HI] were roughly consistent with those from 7r ± 7r ± pairs. 
For hadron-hadron scattering R was usually about one fermi. For nucleus- 
nucleus scattering it was somewhat bigger, perhaps 3 — 4 fm. In both cases 
cr was about one fm, perhaps a bit less. According to most papers, though 
not to all, cr < R and the interaction region had the shape of a pancake 
rather than that of a cigar. Within errors, no systematic dependence of 
R and/or r on energy and/or the kind of hadrons involved was seen. The 
collision energies ranged up to yfs = 52.3 GeV [H], which would have been 
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impossible to study using the old parameters 7 which hardly depend on the 
charges at so high energy per particle [oT] . 

The KP parameterization was not the only one used. In particular Biswas 
et al. in ref. jUj found 

C(p 1; p 2 ) = 1 + (0.80 ± 0.10)e- (11 ' 2±2 - 4)Q2 . (51) 

According to this result i?(p, p) = 0.8 ±0.1. This was surprising, as ac- 
cording to both GGLP and KP, as well as to a more general argument mod- 
elled on quantum optics [SI], i?(p, p) = 1. The experimental evidence for 
R{Pi p) < 1 was n °t yet compelling. Moreover, experimental factors like 
particle misidentification and finite momentum resolution could also reduce 
R(p,p). Nevertheless, the problem attracted attention. 

4.4 Further theoretical ideas 

One way to obtain R(p, p) < 1 is to reject the assumption that the pion 
sources associated with space-time points are incoherent. When all the 
sources act coherently, they can be replaced by one common source. Then 
there is no need to symmetrize and consequently R(p, p) = - there is no 
attraction between the momenta of identical bosons. Fowler and Weiner 
[80J, [HI] pointed out that for a partially coherent source one can obtain 
< R(p,p) < 1. A surprising conclusion was that in order to explain the 
small reduction observed by [H] one needs about 50% of coherence. In a 
more quantitative model proposed in [32] the coherent sources were respon- 
sible for 70% of the pion production. A related strategy was to use a more 
sophisticated model for multiple particle production, where the symmetry 
with respect to exchanges of identical pions is built in from the outset. Work 
along such lines was done by Giovannini and Veneziano [HI)], but unfortu- 
nately it gave the prediction R(p, p) = for multiple particle production in 
e + e~ annihilations, which soon got disproved by experiment. 

Implications of the fact that some of the interfering pions were decay 
products of resonances have been considered since the beginning of the KP 
series of papers [H3]> |108j . In the mid seventies, however, it was established 
experimentally that for multiple particle production in hadron-hadron scat- 
tering the fraction of prompt, i.e produced directly and not from decaying 
resonances, pions in multiple particle production processes is small. Typical 
estimates were between 10% and 20% [93]. This suggested that interference 
effects could be used to study resonances, as foreseen earlier |108j . but also 
that studies of the Bose- Einstein correlations, which did not take resonances 
into account, might be unrealistic jOH], |14fij . In particular Grassberger [9.3] 
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produced a realistic model including resonances and worked out many of its 
predictions. He found that the GGLP and KP models require important 
revisions. Since there is a time lag between the production of the prompt 
pions and the production of the pions from resonance decays, the interfer- 
ence between pions of these two origins should produce narrow peaks which 
should be added to the more smooth functions R(q 2 ) obtained neglecting 
resonances. This gives -R(p, p) > 1, but the peaks corresponding to the long 
lived resonances, in particular to the rj and oo, are so narrow that they are 
missed under the experimental conditions and one expects that effectively 
in practice R(p, p) < 1. Another argument in favor of the narrow peaks is 
that the resonances fly away before decaying which makes the region where 
the pions are created larger. It should be kept in mind, however, that if the 
resonances did not move at all, the narrow peaks would still be there be- 
cause of the time lag. Since the resonance fractions depend on momenta, the 
correlation functions C(p 1 ,p 2 ), should depend not only on the momentum 
difference q, but also on the sum of momenta i.e. on K. The pions from res- 
onance decays, especially from the long-lived resonances, are produced in a 
region, where the density of particles is low, which makes doubtful the claim 
that pion production happens only on the surface of the production region. 
Finally, from the explicit formulae it is seen that the correlation functions 
depend on the orientation of the difference of momenta q. For q parallel to 
the beam axis the peak in Q 2 is narrower than that for q perpendicular to the 
beam. This means that in this model the production region is cigar-shaped. 
A proposal jHSl with a great future was to choose for study, besides the longi- 
tudinal (with respect to the beam) component of q, its transverse component 
parallel to the transverse momentum of the pion pair, later called the out 
component, and the transverse component perpendicular to the transverse 
momentum of the pion pair, later called the side component. In the examples 
worked out in the paper R S ide is always the smallest. R out for p = must 
be equal R S ide from symmetry, but with increasing |p T | it increases and is 
already larger than Ri ong at |p T | = 300 MeV/c. Ri ong decreases significantly 
when |p T | increases from 100 MeV to 300 MeV. At the time the experimental 
data were not good enough to check these predictions. The differences be- 
tween the new picture and the GGLP-KP one were so large that Grassberger 
concluded: "If they [the peaks] are not found, this would have serious impli- 
cations for resonance cross-sections. If they are found, however, the currently 
observed 7t~tt~ correlations cannot be the particle physics equivalent of the 
Hanbury Brown-Twiss effect". 

Another influential contribution was the paper by Yano and Koonin |147j . 
They were interested in heavy ion collisions and their starting point was a 
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classical single particle distribution 



D(x,p) 



E da 



(ttVoT 



-1 exp [-a(x ■ P) 2 + bx 2 ] 



(52) 



a d 3 p 

where P is the total initial four-momentum. The physical origin of this 
formula is simple. The measured momentum distribution is multiplied by 
a guessed distribution in space-time and the product is interpreted as the 
phase-space distribution. Since the source function D(x,p) depends on both 
coordinates and momenta it is, of course, unsuitable for substitution into the 
Klein-Gordon equation. The parameters a and b can be rewritten as 



1 (r 2 + r- 2 ) 



-2 



In the overall cms frame P = {y/s, 0) and formula ()52j) becomes 

E da 



D(x,p) 



a d 3 p 
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The corresponding two-body distribution is 



(53) 



(54) 



EiE 2 da c 
a d 3 pid 3 p 2 



d i x 1 d 4 x 2 D(x 1 ,p 1 )D(x 2 ,p 2 )\(j) s PiiP2 (x 1 ,X2)\ : 



(55) 



where 4> p )P (^1,^2) is the symmetrized two-body wave function. In the 
GGLP and KP models symmetrized plane waves had been used . Yano 
and Koonin considered also Coulomb distorted waves and Coulomb distorted 
waves including strong tt~tt~ interactions in the I = 2, L = 0, 2 states. The 
Coulomb and strong corrections turned out to be marginal (below 0.1% in 
the interesting region of q), so finally only the plane waves were used. These 
give 



^(Pl>P2) = ex P 

and in the overall cms system 

fl(Pi>Pa) = ex P 



a 2 2/ r>\2 1 1 2 

r° T {q - p) + 2b q 



1 2 2 ^-22 

- 2 q r o 



(56) 



(57) 



This formula, or formula f|54j) . gives the physical interpretation of the param- 
eters r and r . As compared to the KP formalism it suggests a change from 
the parameter qx to q 2 . Note also the change, with respect to the GGLP 
formula ()21jl . of the sign in the term proportional to q 2 in the exponent. 
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Let us consider now the quantum-mechanical interpretation of the YK 
formulae. Using the plane wave approximation for function cf) s p p2 and the 
fact that function D(x,p) is real, we get by comparison with ([T8j) 

|p( Pl ,p 2 )| 2 = Re J d 4 x 1 D(x 1 ,p 1 )e^ 1 J d A x 2 D{x 2l p 2 )e~^ . (58) 

This cannot be generally true, since the right hand side does not even have 
to be positive. Following YK we assume further, however, that 

D(x,p) = \A(p)\ 2 p(x), (59) 
where p(x) is real non-negative. Then the Re is unnecessary and one can put 

p(p 1 ,p 2 )=A( Pl K(p 2 )< e ^>, (60) 

which formally coincides with (f4*2*|) and thus, is acceptable as a formula con- 
sistent with quantum mechanics. The factorization of the D(x,p) function 
assumed in (joTJj) means that there are no p — x correlations. This is a rea- 
sonable assumption |101j for the low energy heavy ion collisions considered 
in |147j . but not in general. 

When there is no factorization, the YK integrands contain the non- 
factorizable analogues of p(xi)|y4( Pl )| 2 and p(x2)\A(p 2 )\ 2 respectively, while 
in quantum mechanics the integrands contain the non fatorizable analogue 
of p(x)A(p 1 )A*(p 2 ) and its complex conjugate. This can give similar results 
only under very special assumptions (cf. 13QJ ) . In order to see how the YK 
approach can break down when there are momentum-position correlations, 
let us consider the following example 4 . Replace in formula (JH2|) the normal- 
ized Gaussian by 5 A (x — Xp), where A is a constant. Then the integrals in 
formula (1581) can be done and one finds 



|P(P1, P 2 )| 2 = P(Pl> Pl)p(P2> P2) COS {^Q 2 ) > ( 61 ) 

which is evidently wrong since the right-hand side is not positive definite. 
For some more realistic models where the YK method fails because of the 
x - p correlations and for its derivation using the "smoothness assumption" 
see |T3E| . 

4.5 The GKW paper 

A detailed discussion of particle correlations in heavy ion collisions was given 
by Gyulassy, Kauffmann and Wilson (further quoted GKW) [HE] . This is the 

4 This is a simplified version of the models described in refs |127) . See also . 
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second most quoted paper (after GGLP [88] ) on Bose-Einstein correlations in 
multiple particle production processes. The starting point of GKW was the 
Klein- Gordon equation with a classical current on the right-hand side (J27j) . 
but with A interpreted as a Heisenberg pion field. The normalized solution, 
which had been well known, is the coherent state 



\A) = e-^exp ( i J d i P J(p)a\p)j |0>, (62) 
where a^(p) is the creation operator for a pion with momentum p, 

/gii?pt-jp-x 

^ (2^(2^)3)1/2 ( 63 ) 



n= / d 3 p\J(p)\ 2 . (64) 



and 



The use of coherent states to describe Bose-Einstein correlations is standard 
in optics and had been applied to multiple hadron production processes [S] , 
but GKW pushed the analysis much further than their predecessors. 

Let us denote by | J) the coherent state corresponding to the current J(x). 
Since as is well known 

a\p)\J)= iJ(p)\J), (65) 

one easily checks that [HE] in state | J) the multiplicity distribution is Poisso- 
nian with average multiplicity n. For a given multiplicity m the momentum 
distribution is 

m 

N( Pl ,...,pJ = J] Wl 2 (66) 

3=1 

and consequently -R(p, p) = - there is no GGLP effect. Now the authors 
introduce N source currents related by space-time translations. The density 
matrix is assumed to be diagonal in iV and the probability distribution P(N) 
is part of the theoretical input. Thus the overall source current for given N 
is 

1 N 

7 = 1 
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or equivalently 



i N 

j(p) = Up)^J2 eWkt ~ tp ' x > (® 



N 

where J,r(p) is the Fourier transform of J w (x). The factor l/y/N was not 
used by GKW. It makes, however, the limits N — > oo easier to see (cf. e.g. 
[34] ). The probability distribution in space-time for the points Xj is given by 
some function p(x). It is convenient to define its Fourier transform 

pip) = p*{-p) = J d A xp(x)e tpx . (69) 

Note the crucial difference with respect to the KP approach: it is not assumed 
that the constituent sources J(x — Xj) are incoherent. Thus the solution (|6*2*|) 
can be used with the new current. 

Now the single-particle momentum distribution is 

n(p) = \Up)\ 2 (^ + (^}\p(p)\ 2 ), (70) 

where (• • • ) means averaging over the probability distribution P(N). One of 
the important points of the GKW paper is the observation that since realistic 
functions p(p) decrease rapidly with the increase of any of the components 
of p, they used as a guide formula IjHljl . and since E p > m n the second term 
on the right hand side can be neglected. Another argument for the same 
conclusion is that function p(p) is small, when |p| exceeds 1/-R, where R is 
the radius of the region where the pions are produced, while typical values of 
|p| are of the order of m f or more. Therefore, in the large source limit, which 
is a reasonable picture of heavy ion collisions, the second term is negligible. 
It is a correction term which is calculable, can and should be studied, but as 
a reasonable approximation one can omit it. Then the result is that on the 
average N sources contribute just A" times as much as one source, as if the 
sources were incoherent. In the two-body distribution there are six kinds of 
terms, but after simplifying as in the previous case only two kinds survive: 
those which do not contain p(. . .) and those which contain p(q), where all 
the components of the argument can be made small. One gets 

R(PuP2) = (^}\ P (q)\ 2 . (71) 

For large values of (N) the first factor on the right-hand side is approximately 
equal one and the GGLP result is recovered. Since similar results hold for 
any m-pion distribution, one can assume, as a mathematical shortcut, that 
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the constituent sources are incoherent. GKW called this approximation the 
chaotic field limit and pointed out that this limit, though not the finite N 
corrections to it, had been found from essentially the same model in ref. [3lj . 
Because of their careful handling of N, however, GKW were able to derive 
the incoherence of the sources as a reasonable approximation, while in ref. 

it was introduced as an assumption. GKW realized that the interaction 
region depends on the impact parameter of the collision, which they consider 
to be well-defined in a collision. Thus, they recommended either to use 
data corresponding to a narrow range of impact parameters, or to average 
incoherently over the impact parameters. 

GKW estimated that resonance production is not very important in heavy 
ion collisions. Therefore, in order to explain the fact that -R(p, p) 7^ 1 they 
introduced partial coherence jHU]. Generalizing the proposal of Fowler and 
Weiner they introduced the degree of coherence for momentum p by the 
formula 

D(P) = ( "° (P) , r , (72) 

where no refers to the pions produced coherently and n c h to the pions pro- 
duced incoherently (chaotically). This function should be obtained from the 
relation 

R(p,p) = l-D 2 (p). (73) 

The degree of coherence must be known before conclusions about the geom- 
etry of the source can be extracted from the data, because 5 e.g. 



C( Pl ,p 2 ) = l + (l-/J( Pl ))(l-D(p 2 ))p 2 (g) 

+2 v //J( Pl ) J D(p 2 )(l - D( Pl ))(l - D(p 2 ))p(q). (74) 

Only for D = one recovers the result pSjl . The origin of the coherent 
component was ascribed to some collective effects in the colliding nuclei. 
Thus it would be specific to heavy ion collisions. 

The most influential part of the GKW paper was their analysis of the 
final state interactions, i.e. of the interactions taking place after the pions 
had been produced. These interactions can affect the pion momentum dis- 
tribution. Final state interactions were considered before [§3], |147j . but 
always with the conclusion that they are unimportant. GKW proposed a 

5 Assuming that pip), Jo(p) and Jchip) are real. 



27 



model for these interactions and when it later became used by experimental- 
ists the corrections were often found significant. There are two kinds of final 
state interactions. The pions move in a single particle potential including 
the Coulomb potential from the positive charge of the nuclear cores. More- 
over, the pions interact with each other. Here one can consider the Coulomb 
interactions between the charged pions and their strong interactions |147j . 
Taking into account the single particle potential, say V(x), is comparatively 
easy, because it does not affect the independent particle approximation. One 
assumes that the potential V(x) does not produce particle pairs and does 
not support bound states. Then it is enough to replace everywhere the pion 
plane waves, which are solutions to the free particle Klein-Gordon equation, 
by the corresponding solutions of the Klein-Gordon equation with the po- 
tential V(x). We will denote these solutions ip p (x), where the subscript p 
means that this function corresponds to the plane wave describing a free pion 
with momentum p, i.e. that for x 2 — ■> oo it tends to this plane wave. With 
the additional assumption that the range of the single particle potential V 
is much larger than m" 1 , one finds the single particle density matrix 



From this density matrix the correlation function is built as usual, using for- 
mula (fTHj) to construct the necessary elements of the two-body pseudo density 
matrix. The density matrix may look similar to the density matrix from the 
KP model with extended sources (JHEJ), but the physical interpretation is very 
different. In the KP model the functions (x\i/j(xq)) can be normalized to one 
and describe the shape of a source with a finite space extension. Here the 
functions if) p (x) are almost plane waves extending over all space. They are 
only slightly deformed in a finite region of space due to the potential V(x). 
The effect of coherence is calculated as before, but using the distorted func- 
tions pv (P)P')j ^ov(p) and n c /iv(p)- The distorted functions are calculated 
with the plane waves replaced by the functions tp p (x). In practice this part 
of the final state interaction does not affect much the distribution of the mo- 
mentum difference q, because both bosons are pushed in the same way by 
the potential. 

The two-pion interactions spoil the independent particle approximation 
and, therefore, are much more difficult to include. GKW considered various 
approximations. They pointed out that the distorted wave Born approxima- 
tion would be the best tool, but it is rather complicated. The main effect is 
to change the wave function of the pion pair from a symmetrized product of 
plane waves to a two-body Coulomb wave function. Using the nonrelativistic 
Schrodinger equation, factoring out the motion of the pair as a whole and 




(75) 
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neglecting in the relative coordinate space the changes of the wave function 
in a region of radius O(R) where R is the radius of the interaction region 
around the point xi — x 2 = 0, one finds (HEj that the pseudo density matrix 
element p(p 1; p 2 ; p 1; p 2 ) gets multiplied by the "well-known Gamov factor" 

nt \ 2lT71 am * (tr\ 

G (Pi'P 2 ) = ^T3T' ^IpT^' (76) 

where a is the fine structure constant. According to the theory the measured 
two-particle distribution contains this factor. Thus, the suggestion is that 
experimentalists could correct their results for Coulomb interactions by di- 
viding the measured functions p(p 1; p 2 ; Pi, P 2 ) by the corresponding Gamov 
factor. 

Let us look at the numbers. The Bohr radius for a pair of identical 
particles of mass m and unit charge is 

a B = — , (77) 

am 

where a = 1/137. This is 1.96 MeV™ 1 or 386 fm for charged pions and 0.56 
MeV -1 or 110 fm for charged kaons. Since the Bohr radii are much larger 
than the interactions radii R, which do not exceed several fermis, it had 
long been believed that the Coulomb final state interactions can be safely 
neglected. From the formula for the Gamov factor it is seen, however, that 
the dimensionless parameter is not R/a B , but l/(qa B ), which can be large 
when q is small. Moreover 

Air c , 

27T77 = = -. 78 

qa B q 

The coefficient c is 25 MeV for pions and 7 MeV for kaons. It is enhanced 
by the Air factor. 



5 Refinements and criticism in the eighties 

5.1 New experimental results and problems 

In the eighties the Bose-Einstein correlations were studied in further pro- 
cesses. They were found in \ip |2Zj and vD jTHj collisions and studied in 
detail in e + e~ annihilations jl], [Hj, [T3], j3U], |104j and 77 scattering |104j . 
The radii of the hadron production regions in these processes, as compared 
to hadronic interactions, were found somewhat smaller, as expected. The 
values of R(p,p), however, were not reduced. This contradicted the predic- 
tion from some models (cf. e.g. [BE]) that e + e~ annihilations should be much 



29 



more coherent than hadronic interactions and consequently exhibit little or 
no GGLP effect. In fact, it made the interpretation of A as the coherence pa- 
rameter doubtful. Correlations between charged kaons have been found and 
studied jSj. In general, there was a spectacular improvement in the quality 
of the data. 

As the experimental data improved the uncertainties and ambiguities in 
the analysis of Bose-Einstein correlations became apparent. In ref. jTI] the 
important result jUj that R(p,p) ^ 1 was confirmed. The authors used a 
formula similar to the KP formula (j47jl 



where I(x) = 2J\(x)/x. The introduction of the factor A into the interfer- 
ence term significantly improves the fits to experimental data and was soon 
followed by most groups. Note that in such parameterizations -R(p, p) = A. 
Where we write in the following R(p,p), many writers would have writ- 
ten A. The interpretation of this coefficient, however, remained doubtful. 
Experimental uncertainties, coherence and resonance production were all in- 
voked. The left-hand side is here the ratio of the distribution for pairs of like 
sign pions to the corresponding background distribution of pion pairs. The 
strongest result in ref. |7Tj was that for pairs of like sign pions produced in 
7i + p interactions at 16 GeV/c incident momentum: A = 0.49 ± 0.02. Sta- 
tistically this is about 25 standard deviations from one, but changing the 
background from pairs of unlike charge pions to another background, which 
looked at least as reasonable though later it was criticized [T], the authors 
obtained A = 0.88 ±0.02. Introducing the Gamov factor, which was not done 
in ref. jTI], further significantly increases the fitted values of the parameter 
A |SSj, [HI]- |148j . Incidentally, introducing the Gamov factor makes the fits 
worse. In a particularly bad case |148j the confidence level dropped from 
65% to 0.1%. This explains why not all the groups were willing to include 
this correction factor. The uncertainty in A was particularly striking for pp 
annihilations at rest [TTj- With the background of unlike sign pion pairs, 
A = 1.20 ± 0.08 was the biggest among the three A-s for the reactions con- 
sidered in the paper. With the other background it became 0.63 ± 0.05 - the 
smallest of the three. For Q 2 decreasing towards zero function R(pi,p 2 ) was 
found to increase faster than the Gaussian fit suggested. This was described 
either, following |121j . by including in the fit another Gaussian with a larger 
value of R [H], [2], [25] . |3*Ij . |116j . or by using exponentials in Q instead 
of Gaussians [3], |116j . |123j . Taking into account this additional rise again 
increased the value of A. The uncertainty in A is seen to be large and mainly 




(79) 
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systematic. It was not clear how to reduce it or even how to estimate it. This 
has made progress in the understanding of the origin of this factor difficult. 
The key problem is a good choice of the background distribution. 

Let us note some new ideas concerning the evaluation of the background 
distributions. In order to eliminate the effects of resonances and their re- 
flections, the measured distributions for like meson pairs and those for the 
background were divided by the corresponding Monte-Carlo distributions 



(TU, [IEJ, gjj. A factor (c + c 2 Q 2 ) 0, jUJ , |2Z|, |HH|, [HI!, (co + c x Q) @j 



(with constant Cj-s), or similar [TH], [S], were introduced into Co(pi,P2) to 
keep the correlation function close to one in the region, where Bose-Einstein 
correlations are no more expected. The groups using mixed backgrounds re- 
alized that, in order to get a good mixed background, members of the pairs 
of events supplying the background pairs must be similar. Depending on the 
process this could mean similar multiplicity (cf. e.g. |T|, (SHj, jSHj), rotating 
one of the events to make the sphericity axes coincide (cf. e.g. j3j) etc. The 
mixed distribution was sometimes 0], |148j improved as follows. From (|%8Jl 



where M' is a suitable normalizing factor. In first approximation R(pi, p 2 ) = 
and this is just the single particle distribution as used to build the ordinary 
mixed background. One can do better, however, by substituting into the 
integrand the function -R(p!, p 2 ) calculated in the first approximation. After 
a few iterations one gets a single particle distribution which satisfies the 
equation above. This, however, is a small correction. Note that the product 
of such single particle distributions still does not contain the two-body non 
Bose-Einstein correlations between identical particles, which should be there. 

Since the results for the life-time r were also usually found with large 
uncertainties, most of the discussion concerned the radii R of the interaction 
regions. As a crude characteristics one could say that within the large un- 
certainties all these radii are similar except for the much bigger radii found 
in heavy ion collisions. Nevertheless, it was possible to find some interest- 
ing regularities. Indications were found jS], [HI] that when identical kaons 
are used, the radii come out smaller than the corresponding radii obtained 
when identical pions are used. This could mean a decrease of the effective 
(measured) radius with increasing mass of the identical particles. 

The increase of R with increasing particle multiplicity or particle density 
in pseudorapidity was at first somewhat confusing. For hadron-hadron colli- 
sions, when increasing the particle multiplicity by selection or by rising the 
energy, an increase of R was visible, but only at energies yfs > 30GeV [5], 



one finds 
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0, jZj, [TT] . For e + e~ collisions at energies >/i < 30 GeV no increase 
was seen. For heavy ion collisions increasing the mass number A, or reducing 
the collision parameter gave an increase of R even for energies of the order 
of 1 GeV [El,|in|, |124j . What counts here is the mass number of the 
projectile, provided the target is heavy enough 0. This has been nicely il- 
lustrated in ref. [HHl, where for pp and p + Xe scattering the radii are the 
same within errors, while the particle density in rapidity changes by a factor 
of two. Refs. jH2j, |102| reported for 16 + Au collisions a strong increase 
of Rt (alomst by a factor of two) for rapidities of the pion pairs close to the 
cms rapidity, but this was not confirmed. From a thorough compilation of 
data for heavy ion collisions it was concluded [35] that the measured radius 

R increases linearly with the cubic root of the mass number of the lighter of 

i 

the two colliding nuclei. A fit gave R = 1.21Ap. Some evidence for the de- 
crease of the measured radius with increasing momentum [23] and transverse 
momentum j^U] of the pion pair was also reported. 

Another interesting observation was that for e + e~ [TSJ, [3U], P£)4J and 
hadron-hadron [Zj, [Hj, [SB] interactions, R(pi,p2) could be well approximated 
by a function of Q 2 only - in first approximation by a Gaussian as proposed 
by GGLP. This has in particular two important implications. Contrary to 
the Ansatz made for heavy ion collisions in ref. |147j . for e + e~ annihilations 
the enhancement due to Bose-Einstein correlations is large for Q 2 small even 
if g 2 an d q 2 taken separately are large. The other conclusion is that the 
distribution of q in the rest frame of the pair is spherically symmetric. In 
the KP formula the increase of either q 2 or q\ reduces R{pxiV%) like in |147j . 
but because of the absence of the dependence on q\\ (of the component of q 
along the direction of p : + p 2 ) the fits are much better than in the YK case. 
The reason is that Il5j,[ll 



e * q \\ + H^TT) q * w 1, (81) 

where 7 is the Lorentz factor for the transformation from the cms to the rest 
frame of the pion pair. Introducing this factor into the KP formula for the 
correlation function is approximately equivalent to a change from q^ to q 2 
accompanied by the introduction of a term, which increases with increasing 
go- Thus the KP formula can be rewritten so that it becomes similar to the 
GGLP formula. 



5.2 String model 

For e + e~ annihilations, where the final state is believed to evolve from a sin- 
gle string, the effects of Bose-Einstein correlations were calculated [SI] , JZ| 
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in the framework of a slightly extended [21] string model. The string 
model is closely related to the realistic and highly successful LUND model 
|16j . Therefore, its predictions concerning the distribution and properties of 
particle sources must be taken seriously. Already the qualitative results were 
striking. There are strong correlations between the position of the source and 
the momentum spectrum of the particles it emits. As a result, the volume 
determined from the study of Bose-Einstein correlations is not the volume 
of the interaction region. Only particles with similar momenta contribute to 
the GGLP effect. Therefore, the measured volume is the volume of the re- 
gion, where pions with similar momenta are being produced and not the total 
interaction volume. This has little effect on the transverse dimensions, but 
the measured longitudinal dimension becomes almost energy independent, 
while the true length of the interaction region is believed to grow roughly 
like a/s. In order to introduce x — p correlations into the KP sources |113| it 
is enough to remove the assumption that the sources are equivalent [Slj, but 
to find a quantitative relation between the string model and the KP model 
generalized in this way is not easy. The attempt in ref. [SSj does not look 
very convincing. Thus the geometrical interpretation of the measured vol- 
ume is somewhat doubtful. The one-dimensional string model, which is more 
constrained than the three-dimensional version, suggests that the space-time 
distribution of sources should depend mainly on the variable r 2 = t 2 — z 2 . 
A natural extension to three dimension is to assume that this distribution 
depends on r 2 = t 2 — r 2 , which according to the standard method of analy- 
sis corresponds to a correlation function depending only on the variable Q 2 
[51 J . This is supported by experiment. The model suggests also that the 
distribution in Q is more peaked at small Q and has a longer tail at large 
Q than a Gaussian. At the quantitative level the comparison of the model 
with experiment is difficult and requires the use of Monte Carlo programs 
and additional assumptions e.g. about resonance production. Nevertheless, 
good agreement with the data from ref. jlj was obtained [Hj, [S2], though 
after an interesting detour. At first it seemed [U| that the introduction of 
long lived resonances ruins the agreement. Then it was pointed out that 
the amount of rj (0.41 per event) commonly assumed at that time was only 
a guess and probably a very bad one. Neglecting rf production good agree- 
ment with experiment was obtained. Later experiment confirmed that the 
production of 77' in e + e~ annihilations is indeed small. 
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5.3 Wigner functions, emission functions and covari- 
ant currents 

It has always been a problem how to put together the available informa- 
tion about the space-time structure of the interaction region and about the 
momentum distribution of the produced particles. Usually in quantum me- 
chanics one works in the momentum representation or in the coordinate rep- 
resentation, but not in both simultaneously. KP suggested to consider the 
space-time distribution of the points labelling the incoherent sources and 
a wave function in the momentum representation, the same for each source 
|113j . Yano and Koonin |147j introduced a distribution D(x,p) without both- 
ering about its interpretation in quantum mechanics. Pratt |135j proposed 
to use for the sources the Wigner function (cf. [Iflflj and references quoted 
there). The Wigner function W(X\, K x ) can be defined by either of the two 
equivalent formulae: 

W(Xi,Ki) = y^p^i + qiAKx-qiAtOe^ 1 (82) 

/d 3 Y 
t^4p(Xi + YxAXx - Y 1 /2,t 1 )e- iYl ' Kl . (83) 

Inverting the first equation one finds 

p(Pi, Pi, *i) = j rf 3 Xe-^ x W(X 1 , Kx). (84) 

After freeze-out, the time evolution is by definition that for free particles. 
Therefore, it is convenient to factor out from the density matrix the free par- 
ticle evolution factor [136J, or equivalently to go into the interaction picture. 
This yields 

Pi(Pi,p'i,t) = J ^Xe^WHX^Kx), (85) 

where in the exponent stands the Lorentz invariant product of four-vectors. 
Further the subscripts I will be skipped. 

Now come the characteristic assumptions of this approach: the sources 
are labelled by their freeze out time to, and the summation over the sources, 
assumed incoherent, is reduced to an integration, with a suitable weight 
factor, over the freeze out time t . Since for any source p/(p, p',t) does not 
depend on the time t for t > to, where to is the freeze out time of this source, 
one can put t = to on the right hand side. Then the formula for the full 
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single particle density matrix at any time t larger than the latest freeze out 
time, reads 

p(Pi,p'i) = J d'Xe^SiX,^), (86) 

where the matrix element on the left hand side is time independent and the 
time component X on the right hand side is understood as the freeze-out 
time to of the source labelled to- As seen from this derivation function S, 
which will be further referred to as the emission function 6 , can be understood 
as a product of the Wigner function for the particles produced by the source 
to and of the density factor necessary to convert the summation over the 
sources into the integration over t . The point is, however, that besides this 
emission function there is an infinity of other emission functions which, when 
substituted into this formula, yield the same density matrix and therefore are 
equally acceptable. In particular a small source with a large life-time can give 
the same correlation function as a large source with a short life-time [136 . 
The statement |135j that S(x, p) "can be identified as the probability of 
emitting a pion of momentum p from space-time point x" suggests a way of 
thinking about the source function which has become very popular. What 
can be proved, however, is much less: the Wigner function W(x, p) when 
integrated over space (d 3 x) gives the momentum distribution at time t and 
when integrated over momenta (d 3 p) gives the distribution in space-time 
[100 . Nevertheless, it soon became customary to discuss the source function 
using Pratt's picture. 

When sources freezing out at different times contribute coherently the 
relation of the emission function to Wigner functions becomes rather tenuous 
|151j . In this case it seems better to relate the emission function to the 
history of the sources rather than to some states of the system |139j .[52]. Let 
us define the emission function in space-time by the formula (c.f. formula 

m) 

S(X 1 ,Y 1 ) = Y,(Mxi)J!(x' l )). (87) 

i 

Then according to Shuryak's formula (|3*3*jl 

p(Pi, pi) = j S(X, Y)e iKlY+i « lX d 4 Xd% (88) 

6 Also known as Wigner function, pseudo Wigner function, source function etc. The 
name emission function was used in ref. |lHf>| . 
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Formula (JSfij) is reproduced for 



S(X 1 ,K 1 ) = j S(Xi,Y)e iKlY d 4 Y. (89) 

This is an unambiguous definition of the emission function. The emission 
function is real, because complex conjugation is equivalent to a change of the 
sign of Y and this can be compensated by a change of the integration variables 
from Y to —Y. Note, however, that in order to calculate this emission 
function it is not enough to know the state of the system at some time, or 
even to know the state of the system at every time. Phase relations among 
the states generated at different times are also important. The quantity being 
averaged depends on two times, and not on one as would be case for state 
vectors, density matrices or Wigner functions. 

Using the formula for the single particle density matrix and formula ()18j) 
one finds 

1/ d A XS(X,K)e iqx ? 

c(Pij P2) = 1 + fMs(x, Pl )fM's(x', P2 r (90) 

Often it is considered an acceptable approximation to replace each of the 
momenta in the denominator by K. Then the formula takes a simpler form 



C( Pl ,p 2 ) = l 



j d*xS{x, K)e iqx 2 



(91) 



/ d*xS(x,K) 

The second term on the right hand side is the function -R(p 1; p 2 )- 

An alternative approach has been worked out by Gyulassy and collabo- 
rators [2fl|, |106j . jHZj, [ESI], [US!!- ^ * s an extension of the labelled wave 
packets method of KP [US]. Suppose that each source (wave packet) is 
labelled by a space-time point x s and a momentum p s . The energy corre- 
sponding to this momentum is obtained from the condition p 2 = m 2 . Then 
one assumes that 



p(Pi,P 2 ) = <e i9 > Ps (Pi,P2)>, (92) 

where the averaging is a summation over the sources s. The new point, as 
compared to the KP formula (|42|). is that the single source density matrix 
on the right-hand side depends on p s . A particularly important special case 
is when 

PpAPi,P2) = J{—— )J (— -)• ( 93 ) 
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Interpreting p s as the average momentum of the pions emitted by the source 
this means that each source looks the same in its rest frame, while for KP 
each source looked the same in the overall frame of the event (except for a 
space-time translation). Also the condition p 2 = m 2 means that the labels 
correspond to the emitted pions rather than to a heavy source. Further, 
it is assumed that the averaging over the sources can be expressed as an 
integration with weight function D(x s ,p s ). Thus 



This version of the model is known as the covariant current formalism. A 
more general model, where the single source density matrix does not factorize, 
i.e. does not correspond to a pure state, and the second equality does not 
hold, was considered in |13()| . The authors proposed 



where N is the normalization constant, while Ax 2 and Ap 2 are constants, 
variances of x and p respectively, constrained by the Heisenberg uncertainty 
relation Ax 2 Ap 2 > 1/4. This density matrix factorizes if and only if the 
terms proportional to pip 2 in the exponent cancel, or equivalently when the 
product Ax 2 Ap 2 has the smallest value allowed by the Heisenberg uncertainty 
principle, i.e. for minimal wave packets |130j . The factors are then functions 
of (p s — Pi) 2 = 2m 2 + 2p s pi, i = 1, 2 and the density matrix can be written 
in the form (jHSJ) 

Both approaches described here use simultaneously momenta and coordi- 
nates. Their interpretations, however, are different. In the Wigner function, 
or Shuryak, approach X and K are arithmetical averages of the arguments 
of p(x, x') and p(p,p') respectively. In the wave packets approach they are 
labels of the sources, which may, but do not have to, be the centres of the 
wave packets. Single particle energy is no problem since the particles are on 
mass shell. Time in the wave packet formalism is just one more label. In the 
Wigner function approach its meaning depends on the details of the model. 




(94) 



Pp.(Pl>P2) = NeX P 



{p s -{pi+p 2 )/2) 2 { Pl -p 2 ) 2 Ax 2 
2Ap 2 2 



(95) 
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5.4 Difficulties with the final state interactions 



Final state interactions were further studied. It had been known (cf. Section 
4.5) that the two-body Coulomb interactions are controlled by the dimen- 
sionless parameter £ c = R/as 136 , where as is the Bohr radius and R is the 
radius of the interaction region. For pion pairs a# ~ 386 fm. Consequently 
£ c ~ 0.01 for heavy ion collisions and much less for the other cases. There- 
fore, it had been believed for some time that Coulomb correction cannot be 
important. This conclusion is not justified, because the factor R in the small 
parameter £ c enters as an estimate of l/|q|. In the very small |q| region this 
is a gross underestimate. Moreover, as seen from the Gamov factor, there 
are large numerical coefficients. The leading correction in the small £ c limit 
is given by the Gamov factor |96j . This correction is significant and spoils 
the agreement with experiment in the small Q 2 region (cf. e.g. |85j). 

Pratt |136j has used the GGLP model with a Gaussian distribution of 
sources and with the symmetrized product of plane waves replaced by a 
symmetrized Coulomb wave functions. This method of taking into account 
Coulomb interactions had been used before |107j . |120j . but mainly in the 
study of the Fermi-Dirac correlations among nucleons. Pratt found that 
including the next term of the expansion in powers of r], besides the Gamow 
factor, changes the estimated interaction radius by about 20% for R as small 
as a few fermi. 

In order to describe the two-body strong interaction among identical 
charged pions it is enough to consider the s-wave, 1 = 2 phase shift |144j . 
The observed phase shift corresponds to repulsion, as one would expect for 
an exotic two pion system where no resonances can be formed. The key 
observation [33] is that since the interactions are short range (about 0.2 fm) 
they are negligible for pions produced at a distance of one fermi, or so, from 
each other. Consequently, this correction for strong interactions (not to be 
confused with the corrections for resonance production) is negligible for the 
study of Bose-Einstein correlation in heavy ion collisions. For e + e~ annihi- 
lations, on the other hand, this correlation can be important. According to 
a rough estimate [33] it reduces the factor A by at least 30%. Moreover, in 
agreement with experiment, it produces a minimum of the correlation func- 
tion at Q ~ 0.6 GeV, where the correlation function drops below one. In 
spite of this success the quantitative applicability of this model, where only 
an isolated pion pair is considered, is doubtful [33] . For 7r + 7r~ pairs there is 
an attractive strong interaction due to the 1 = s-wave phase shift. There- 
fore, three-body 7r + 7r~7T + interactions give an effective attraction in the 7r + 7r + 
system which, according to rough estimates, could easily overcompensate the 
repulsion due to the 1 = 2 phase shift [SS] . A more careful analysis including 
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resonance production [56 a suggests that the strong interactions must some- 
how cancel, because otherwise, the fully corrected A would come out negative, 
which does not make sense. Another important remark [HI] is that experi- 
mentally, at low Q 2 , the backgrounds formed from unlike sign pions are not 
very different from the mixed backgrounds, where by construction the strong 
interaction among members of pairs is absent. This supports the point of 
view that final state interactions are less important than naively expected. 
A study of two-body correlations for pairs of unlike sign pions could bring 
here 7 interesting additional information. Bowler's final conclusion for the 
strong final state interactions is [SZj "I think there is good reason to believe 
that in multiple production the effects of final state interactions among all 
pairs largely CANCEL and only the propagation and decay of resonances 
modifies the underlying source structure". An exact analysis of such many 
body strong interaction is, however, "impossibly complicated" j^E]- Bowler 
concludes also |55J that "the chaoticity and range parameters extracted from 
experiment have no simple interpretation and should be regarded as purely 
descriptive" . 



5.5 Early models with x - p correlations 

The improved formalism has been used to build models more realistic than 
the early attempts of GGLP and KP. A common feature of these models is 
that they include x — p correlations. Thus Pratt |135| considered 



S(x, p) = 8(r — R)5(t) exp 



E v — vr ■ p 
-7 f 



(96) 



where the temperature T, velocity v and Lorentz factor 7=1/ yl — v 2 are 
constant. The pions are ejected from a sphere of radius R and tend to 
go radially outwards. Thus, contrary to the KP picture, there are x — p 
correlations. Evaluating i?(p 1 ,p 2 ) for this source function one finds |135j 
that the effective (i.e. measured) radius of the interaction region is: 

R eff (K) = ^(ytanhy)- 1 - (sinhy) 2 ] 1 / 2 ; y = E^. (97) 

The expression in the square bracket equals | for |K| = and decreases 
towards zero with increasing |K|. However, in this model we know that the 
true radius of the interaction region is R for all K. Let us stress the impli- 
cations of this important observation. One should study the Bose-Einstein 
correlations at given |K|. A |K| dependence of the deduced radius of the 



7 And in fact latter brought. 
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interaction range implies the presence of x — p correlations. Wrong assump- 
tions about the x — p correlations lead to wrong determinations of the radius 
of the interaction region. For instance in Pratt's model averaging over |K| 
and interpreting the results as if there were no x — p correlations one obtains 
a radius below a/2/3 of the correct (input) value. 

A more realistic model with x — p correlations, inspired by Bjorken's 
inside-out side cascade, was given in ref. |lUb'j in the framework of the covari- 
ant current approach. The authors chose 



D(x,p) 

■ (PsP, 

J — 

m 



S(t - r coshy)(5(z - r smhy)e *t/ r t^ 



a cxp 



PsP 

2mT 



(98) 
(99) 



where a, T, and r are constants. This corresponds to a distribution of 
sources which is uniform in rapidity and concentrated at the longitudinal 
proper time tq = \Jt 2 — z 2 . Each source in its rest frame produces pions 
according to a "pseudothermal model" which is a modification of the Boltz- 
mann distribution chosen so as to make the subsequent integrations possible 
to perform analytically. The corresponding single particle pseudo density 
matrix is [Ml, [1301 



p( Pl ,p 2 ) = aK (V^e-^; (100) 
m 1T + m 2T 



a 



2T 

1 



iT (miT - m 2T ) 



+ 



4T 2 



+ r m lT m 2 T [cosh(yi - y 2 ) - 1] , (101) 



where rriiT are the transverse masses, the rapidities of the two particles 
and Kq is the modified Bessel function. It is interesting that this model, 
extended by the inclusion of resonances which requires the use of a Monte 
Carlo program, explains the data of the NA35 experiment [3*2*] . which had 
been previously explained [H"] by assuming the presence of a quark-gluon 
plasma. Resonance production is the most important correction, because 
other corrections like a spread in the freeze-out proper time tq or a nonuni- 
form distribution in rapidity tend to cancel |129j . This analysis can be gener- 
alized and made more rigorous |130j . Let us also note attempts to use Monte 
Carlo programs to get more realistic distributions of the sources |101j . [66J. 

A rather different approach to Bose-Einstein correlations in hadron and 
heavy ion collisions was developed by Makhlin, Sinyukov and collaborators 
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(cf. j!25j . [2H], |126j . |143j and references quoted there). Their starting point 
was Landau's hydrodynamic model. According to this model the production 
of hadrons in high energy hadronic interactions proceeds in three stages. 
First the content of the interaction region thermalizes. This happens very 
fast. Then the fluid expands and cools according to the laws of hydrody- 
namics. Whenever the temperature of an element of the fluid drops to a 
critical temperature T = T cr , its content gets converted into hadrons which 
have the thermal equilibrium distribution corresponding to temperature T cr . 
This hadronization of a fluid element is also assumed to be a very rapid pro- 
cess. Since at every space point the temperature is a function of time, the 
space-time points where T(x,t) = T cr form a three-dimensional hypersurface 
S cr in space-time and all the hadrons are produced there. In order to get 
predictions one should in principle assume some initial conditions i.e. the 
distribution at the beginning of the hydrodynamic stage, solve the hydrody- 
namical equations with these initial conditions and find the hypersurface S cr 
as well as the velocity distribution of the fluid on it. This would be a very 
difficult task, but there are simple estimates. E.g. one can define S cr by the 
equation t 2 — z 2 = r 2 , where r is a constant and assume that the fluid velocity 
is parallel to the z axis and at space-time point x, t equal z/t. The authors 
suggest that by selecting events with a given heights of the central plateau 
in rapidity one can approximately fix r. There is also a more complicated 
estimate due to Landau. When transverse expansion is included, parts of 
the hypersurface S cr become time like which complicates the formalism, but 
can be dealt with |143j . The hydrodynamic approach got very popular later, 
though usually it just means introducing collective velocities of the sources. 

5.6 General proposals 

Some results more general than specific models have also been given. For 
heavy ion collisions the plausible assumption that the transverse profile of 
the interaction region is given by the overlap of the two nuclear distributions 
explains easily two experimental facts. For a given pair of ions the average 
transverse radius of the interaction region increases with decreasing im- 
pact parameter [33J and for central collisions Rt is given by the radius of the 
smaller of the two colliding ions jB0]- 

In ref. j!17j the authors, generalizing earlier results from |132j . proved 
that when p(p 1; p 2 ) integrated over p : + p 2 can be written in the form 
p(a 2 q^ + g 2 ), where a is a constant, the distribution of the angle 6 between 
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vector q and the axis of the event z is 



a 2 



4>(cos9) = —j-z - - , . (102) 

v ' 2{a 2 + (1 - a 2 ) cos 2 9fl 2 y ' 

This formula, which follows by simple integration, can be useful to distin- 
guish, cigar shaped (a 2 < 1) from pancake shaped (a 2 > 1) sources. Pratt 
|137| and Bertsch [32] reintroduced for the difference of momenta q the co- 
ordinate system which had been used some ten years before by Grassberger 
jnn|- The notation and terminology is from Bertsch [1*2*] : qi for the longitu- 
dinal component along the event axis, q Q for the outward component parallel 
to Ky and q s for the sideward component perpendicular to both the event 
axis and to K. His formula 8 

R{ Pl ,p 2 ) = e -(Rh 2 s +R 2 Q 2 o+Rhl) (103) 

soon became very popular. The main point of both authors was that, if 
the time interval when the pions are produced is long - this could be due 
to a slowly hadronizing quark- gluon plasma but other possibilities were also 
considered - then one expects R Q 3> R s , because the pions produced late 
behave like pions produced far along the out direction. Let us finally note 
a diagrammatic classification of the various terms occurring in the theory of 
multiparticle Bose-Einstein correlation with coherence included |46| . For an 
attempt to estimate the effect of multibody Bose-Einstein correlations on the 
two particle correlation function by direct calculation of the permanents see 

pun. 

Podgoretsky |133j |134| studied the dependence of the parameters R and 
r on the boosts along the event axis. Let us denote the event axis by z and 
the axes in the out and side directions by x and y respectively. Let us further 
assume that in some reference frame the distribution of sources is given by 
a formula similar to the Yano-Koonin formula (|54jl : 



p(x) = ~ — exp 

An 2 R 2 RT 



x 2 + y 2 z 2 t 2 



(104) 



2R 2 2R 2 2T 2 

Then, in the reference frame moving with respect to this frame with velocity 
j3 (Lorentz factor 7) parallel to the z axis, the distribution of q is [133J 

W ~ 1+exp \-q 2 x R 2 - q 2 y R 2 - 7 2 {g 2 (l - /3u z ) - q x (3u x } 2 R 2 - -f 2 {q z (P - u z ) - q x u x } 2 T 2 

(105) 

8 Actually in the paper the coefficient of q\ reads but this seems to be a misprint. 
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Here u is the velocity of the pion pair K./K and the identity 

go = u • K (106) 

has been used. Yano and Koonin |147j got a similar result, but instead of 
using explicitly the Lorentz transformation they expressed the exponent in 
terms of invariants, see formula t)52j). Podgoretsky noticed that whenever 
there is a frame where all the pion sources are at rest and where their dis- 
tribution is invariant under the change of orientation of the event axis, the 
parameters Rl and r must be extremal in this frame. Thus, it could be pos- 
sible to identify the frame, where all the sources are at rest. Early estimates 
|20j . [21] favored the quark frame, but later the much more detailed analysis 
of [2] demonstrated, that there is no clear minimum of Rl as a function of 
Pp/Pit+ an d consequently no privileged frame, where all the sources are at 
rest. On the other hand the parameterization (jl()5j) turned out to be very 
useful. 

6 Conclusions 

In his review talk entitled The GGLP effect alias Bose-Einstein effect alias 
HBT effect at a conference in Marburg in 1990 [90]] G. Goldhaber said: What 
is clear is that we have been working on this effect for thirty years. What is 
not as clear is that we have come much closer to a precise understanding of 
the effect. Let us discuss this statement from the present point of view. 

There is no doubt that GGLP was a great paper. It correctly identified 
the Bose-Einstein correlations as the reason for the small Q 2 enhancement for 
the like sign pion pairs in pp interactions 9 . Moreover, the authors proposed 
the invariant Gaussian parameterization, noticed the possibility of getting 
from BEC information about the interaction region, pointed out that a joint 
analysis of many experimental distributions should be made and mentioned 
the possible importance of resonances. 

The most important message: the possibility of measuring the interaction 
region, however, passed unnoticed. Only after the work of Kopylov, Pod- 
goretsky, and Cocconi in the early seventies it became gradually accepted 
that BEC are more interesting as a tool to study the interaction region than 
as a peculiarity of the distribution of the relative momentum of identical 
pions. This was certainly a breakthrough. Another breakthrough happened 

9 Amusingly, in the special case of pN interactions at low energy there is now evidence 
|84| . |19j . |2H1 that the GGLP effect is due to resonances more than to BEC between 
directly produced pions 
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on the technical side. The GGLP method of modifying the integrand of 
the integral over phase space was applicable only to low multiplicity exclu- 
sive channels, while the method, pioneered by Kopylov and Podgoretsky, of 
modifying the correlation function could be applied at any (sufficiently high) 
multiplicity and also to inclusive distributions (at sufficiently high energy). 

The relation between the properties of the interaction region and the 
data on interparticle correlations turned out to be much more complicated 
than anticipated. Final state interactions, resonances, experimental reso- 
lution and particle identification are all important factors. Even, however, 
when all these technical problems are solved, important sources of systematic 
error remain. Once it is recognized that the hadrons are not all produced 
instantly and simultaneously, the problem arises: how to reconstruct a four- 
dimensional, space-time distribution of sources from the measured three- 
dimensional distribution of momenta. Moreover, BEC can be only observed 
among particles with similar momenta. Particles with similar momenta, ac- 
cording to many models, are produced in a homogeneity region, which is a 
fraction of the total interaction region. Thus, what one measures is the size 
of the homogeneity region and not that of the total hadronization region. In 
other words, it is not possible to draw conclusions from BEC without pre- 
vious knowledge about the position- (in space-time) momentum correlations 
of the produced hadrons. 

The conclusion is that attempts to get information about the interaction 
region from a model independent analysis of the data seem hopeless. It 
is necessary to start with a model and then, within the model, it may be 
possible to use the data on BEC to get information about the interaction 
region. This information may be crucial, like in the case of hydrodynamic 
models, or of marginal interest like in string models. 

To summarize: the idea to study BEC as a way of learning about the in- 
teraction region has certainly made their study much more interesting. Many 
problems not anticipated by GGLP have been identified and some of them 
have been solved. A reliable and non controversial method of getting quan- 
titative information about the interaction region from BEC is still, however, 
not available, as was strikingly illustrated by the failure of the predictions 
made for RHIC. If this is meant by the precise understanding of the effect 
Goldhaber's position is still defendable. 
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